Recurrence relations for spin foam vertices 
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. . . We study recurrence relations for various Wigner 3nj-symbols and the non-topological lOj-symbol. 

0^ , For the 6j-symbol and the 15j-symbols which correspond to basic amplitudes of 3d and 4d topological 

■ spin foam models, recurrence relations are obtained from the invariance under Pachner moves and 
' can be interpreted as quantizations of the constraints of the underlying classical field theories. 

We also derive recurrences from the action of holonomy operators on spin network functionals, 
^ . making a more precise link between the topological Pachner moves and the classical constraints. 

Q ' Interestingly, our recurrence relations apply to any SU(2) invariant symbol, depending on the cycles 

•"7" ^ of the corresponding spin network graph. Another method is used for non-topological objects such 

as the lOj-symbol and pseudo-isoceles 6j-symbols. The recurrence relations are also interpreted in 
terms of elementary geometric properties. Finally, we discuss the extension of the recurrences to 
take into account boundary states which leads to equations similar to Ward identities for correlation 
functions in the Barrett-Crane model. 

cr. 
^ . 

bf), I. INTRODUCTION: SPINFOAMS AND THE BC VERTEX 

Spin foam amplitudes A provide a non-perturbative and background independent definition of the path integral for 
^ \ general relativity. Considering a four-dimensional manifold with a triangulated boundary S provided with a given 
^'j^ . classical (discrete) metric qab (on E), then the spin foam proposal is: 

(N ■ A[qab]^ K[qab]= I?ff^,e^^°^[9-I. (1) 

■I qab 

In particular, a key property of the right hand side is to project on the kernel of the spatial diffeomorphisms Tia and 

■ Hamiltonian constraint Ti present in the action, i.e. to satisfy the set of formal equations 

> 

X, .. 

, general relativity. Therefore one of the requirements of a consistent spin foam model is to implement some version of 
these equations. In its present formulation, a spin foam model is defined as a sum over bulk triangulations, 

= ^ A^c^g.fa] (3) 

V 

where v is the number of 4-simplices of the triangulation and A a dimcnsionless coupling constant - the notation v 
comes from the fact that 4-simplices are dual to vertices in the 2-complex dual to the triangulation. This justifies the 
name "vertex expansion" for the sum The amplitudes c„ here include a sum over the triangulations with the same 
V and are built with an elementary amplitude assigned to each 4-simplex. The amplitudes are algebraic quantities 
from the representation theory of the local Lie group of general relativity. A key result of this construction is to 
represent the boundary metric in terms of discrete quantities, i.e. piecewise flat metrics with discrete values, typically 
half-integers. This provides a match with the labels of spin network states, the complete basis of the kinematical 
Hilbert space of LQG. Spin foam models can then be used to give transition amplitudes to spin networks. To support 
the conjecture that the amplitudes represent quantum gravity, one can study equations like ^ and in particular link 
the operator H to some version of the Hamiltonian constraint of general relativity. A priori, the connection with LQG 
suggests that 7i should be an operational version of Ashtekar's Hamiltonian Ti. ^ EEF, however such a link is still 
missing for existing spin foam models (see however [H [3, [3, 0, Q ) . 



T-iaiqab, S/Sqab) K[qab] = 0, Hiqab, 5 / Sqab) K[qab] = 0. (2) 
These Whceler-DeWitt equations encode at the quantum level the invariance under space-time diffeomorphisms of 
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On the other hand, there are semiclassical approximations of spin foam models which suggest a different direction to 
interpret Ti. In these approximations, the piecewise fiat metrics qab are seen as Regge geometries, and the amplitudes 
Cii[9ab] are given by exponentials of the Regge action, thus providing a clear link with general relativity. In this 
approximation, ^ should be realized with H being an operatorial version of the Hamiltonian for Regge calculus. 
Unfortunately the construction of the canonical framework for Regge calculus is still under development 0, S, 0, 0] , 
and there is no clear proposal for the Hamiltonian. The difficulty of representing the Hamiltonian constraint in these 
variables has two origins. The first one comes from the discretization of the boundary manifold. Schematically, the 
form of the classical constraint is 

n^^K'^ + R, (4) 

where K and R are respectively the extrinsic and intrinsic curvature of the 3d manifold. In a classical simplicial 
context, these two quantities are discretizcd respectively on n — 1 and n — 2 simplices. Therefore different realization 
of H can be a priori envisaged. The second difficulty comes from the fact that on top of being a Regge geometry 
as opposed to a continuum one, Qab in © is also represented by discrete quantities such as half-integers. Hence the 
version of ^ we are after should be a difference equation. This second difficulty offers a new light on investigating 
the issue. In fact, difference equations arise naturally in the representation theory of Lie groups used to build the 
amplitudes in ([3|), thus it is tempting to look there for an implementation of the idea of representing the Whceler-de 
Witt equations as difference equations. 

There are two instances where this idea is at least partially realized. The first one is loop quantum cosmology 
(see @ and references therein) . There the mini-superspace setting allows a construction of the Ashtekar Hamiltonian 
Ti ~ EEF (which unambiguously lives on n—2 simplices) precisely as a difference equation. In turn, the exponentiated 
action has been recently argued 0] to give rise to a vertex expansion of the form 

The second instance is the Ponzano-Regge model for 3d Euclidean quantum gravity. Recurrence relations for the 
6j-symbol are well-done [1, Q . They can be used to derive the asymptotics of the 3d quantum gravity amplitudes at 
large scales and for fast numerical computations of the 6j-symbol. We begin our paper with a discussion of 

recurrence relations for this model and their connection to Wheeler-de Witt equations. Section |TT1 We discuss various 
techniques to derive recurrence relations, which allow us to reobtain the existing ones and to produce new ones. One 
of these methods is very standard, making use of the topological invariance of the model. A second method gives a 
clear relation between the classical constraints and the 1-4 Pachner move, using the action of holonomy operators on 
spin network functionals (as a simple case of the analysis of [HI). Finally, for the special case of isosceles 6j-symbols, 
we present a new technique which will be generalized to the lOj-symbol of the non-topological Barrett-Crane model. 

Then we tackle the issue of four dimensional models. In Section [Till we consider Ooguri's model This is a 

topological model corresponding to the quantization of BF theory and not to quantum gravity, but is often used as 
a starting point to build spin foam models for 4d gravity (see for the construction of the most recent spin foam 
models). We show that the topological invariance again encodes difference equations for the SU(2) 15j-symbols. Like 
in 3d, such relations provide a connection between the classical flatness constraint F ~ and the Pachner moves. 

For the above topological models, this shows that the invariance under Pachner moves contains the dynamics 
generated by the underlying classical symmetries. These symmetries are such that for a given topology, all bulk 
triangulations contribute to the sum ^ with the same amplitude (upon regularization) . Thus, the projection onto 
the kernel of the constraints is still realized when the sum (jSj is restricted to a single triangulation, say to the first 
order in A. However, this is not the case anymore for non-topological models. In such models, we expect the classical 
symmetries to be restored, and the spin foam amplitude ^[^ab] to project on the kernel of the constraints only through 
the sum over the bulk triangulations for a given boundary. But we do not know at the present day how to get closed 
expressions for non-trivial triangulations. 

Nevertheless, this is not the end of the story, and it is certainly interesting to look at difference equations even 
for the first order in A. If we formally think of the projector A as the exponentiation of the Diffeomorphism and 
Hamiltonian constraints, and of the sum over spin foams as a sum over histories of spin networks, then we expect the 
first order of ([3]) to be related to the matrix elements of the constraints. As already mentioned, there are difficulties 
in representing the constraints. Thus, recurrence relations are certainly useful to probe their basic properties. In this 
regard, we emphasize in this article the geometric interpretation of the recurrence relations we obtain. In particular, 
for the topological models mentioned above, we found that the basic recurrences are generated by gluing a flattened 
simplex to an initial one. From the point of view of the latter, it results in an elementary displacement of a point, 
which induces shifts of lengths or areas. 

Recurrence relations could also put restrictions on the form of the admissible higher orders and help the renormal- 
ization process. Indeed, if a correction of order n violates a difference equation, it implies that it does not fulfill some 
geometric property of the initial amplitude, and thus might be irrelevant for renormalization. 



With these perspectives in mind, we finally focus on the Barrett-Crane model [21| in Section |Vl This is a non- 



topological model which assigns an amplitude, called the lOj-symbol, to any 4-simplex whose areas take discrete values. 
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Although it has been shown not to provide a correct quantization of general relativity, it captures some important 
features provided we only look at a single simplex. We derive a difference equation of order 4 for the lOj-symbol, 
which is solved in the large area limit by the cosine of the Regge action. This leads to a nice geometric interpretation 
showing that the recurrence probes the closure of the 4-simplex. Then, we propose to look at difference equations not 
for A alone, but together with a boundary state. Such a state should be thought of as a physical state, satisfying the 
Wheeler-DcWitt equations for discrete metrics. By peaking such a state on a well-defined boundary geometry, we 
may expect to suppress the non-geometric parts of the recurrence relations, which are for instance associated to the 
degenerate configurations of the lOj-symbol 0, H^. This is precisely the strategy used in semi-classical computations 
of observable correlations in the Barrett-Crane model [1^. In this context, the recurrence relations moreover translate 
into some Ward- like identities for insertions of observables, as was recently explained in the 3d case in [l0| . 



II. RECURRENCE RELATIONS FOR {6j}-SYMBOLS 

It is well-known that the 6j-symbol satisfies the following recurrence relation, 

A+iU) ( . T ^'\+Ao{j) (^.1 ^A+A^,{j) ( . ^1 = 0. (5) 

+ l 32 Jsj ^ ' \jl 32 33} ^ ' 32 33} ^ ^ 

The exact values of the coefficients are not important to our considerations here, but the interested reader can found 
them in the Appendix. This recurrence relation is a second order difference equation for one of the spins (here chosen 
to be ji, but any other can be chosen suitably changing the coefficients). This is one of many recurrence relations 
satisfied by the 6j-symbol (see e.g. for examples), which can all be derived from the Biedenharn-EUiott identity. 
The latter (together with the orthogonality relation and a basic associativity property, see @) fixes uniquely the 
properties of the 6j-symbol, hence it can be seen as its defining relation. In particular, it can be used to prove that a 
state sum model for three dimensional quantum gravity based on a tetrahedral amplitude given by the 6j-symbol is a 
topological invariant. The Biedenharn-EUiott identity turns out to be precisely the statement of the invariance under 
the so-called 2-3 Pachner move for that model. It thus intimately related to the special homeomorphism symmetry 
of three dimensional gravity. Since this is in turn the symmetry imposed by the classical constraints, one can argue 
that the recurrence relations capture a quantum version of the Wheeler-de Witt equation ([2]). 

A connection between ([5]) and a partial differential equation can be made taking the continuum limit. Let us define 
f{ji) — y^l2TTV{£i) {6j}, where V is the volume of the tetrahedron t whose edge lengths are given by £i ~ ji + 1/2, 
and corresponding dihedral angles by Oi{ji). If we recale homogeneously the spins, i.e. ji i— > Nji, and send oo, 
the relation ([5]) is turned into the following second order difference equation Q : 

[A,, + 2 - 2cos0i(jO] , . \ ,. f{3i) = 0, (6) 
^sm6'i(j() 

where the difference operator A is defined through: A/(a;) — f{x + 1) + f{x — 1) — 2f{x), and 9i is the dihedral angle 
between the two triangles meeting at the edge carrying the spin ji in r, computed from the edge lengths In the 
large spin limit this difference equation is approximated by a second order partial differential equation whose solution 
leads to the known asymptotics of the 6j-symbol in terms of the Regge action^ S-s. [1, [IJ , 

31 32 J3| ^ 1 eos + . (7) 

The large spin regime then provides a notion of semiclassical limit in which the spin foam amplitudes are approximated 
by quantum Regge calculus, see psi. l29j. 

The asymptotic expression (O of the recurrence relation suggests to look for a discretization of the classical Hamil- 
tonian constraints on the links of the triangulation, as pi — 9i {ji ) 0, 0] ■ This can be seen as the realization of all the 
constraints (not just the Hamiltonian (|4])) as = 0, a form peculiar to 3d gravity. 

However, apart from this allusive analogy, a more explicit connection between ([5]) and a discrete version of the 
WdW equation has never been achieved so far (see also [lB|)- The difficulty lies in the difficulties with a canonical 
version of regge calculus, and on the question of whether the constraints should be discretized. 



^ In general, the solution is combination of Airy functions, which reduces to (0 when V{ie)'^ > 0. 
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A simplification can be obtained studying homogeneous configurations, i.e. taking all the spins equal to j and 
studying what happens under variations of j. No analytic recurrence relation can be built in terms of j alone (we will 
come back to this point below). Nevertheless in the large spin limit we can use the explicit solution ^ to derive the 
following equation, 

[a| + (60o)'] fU) = 0, (8) 

where 9o = arccos(l/3) is the dihedral angle of an equilateral tetrahedron. A natural interpretation for this equation 
is a cosmological model. The amplitude associated to the tetrahedron can be seen a the Hartle-Hawking state for a 
two-sphere 5^ discretizcd by four triangles. In fact any higher triangulation in the bulk would be divergent because 
of a gauge symmetry, and fixing the corresponding gauge degrees of freedom following [3^ would reduce the bulk 
triangulation to just the trivial one of a single tetrahedron. This is of course a model which is way too simple to be 
meaningful,^ and one should at least extend these considerations to a more complicated boundary triangulation and 
the resulting spin foam amplitude. A systematic analysis of recurrence relations could be useful to this end. 

We now present a new type of recurrence relations, that can be obtained using an integral representation of the 6j 
-symbol. The latter can be given for a special "isosceles" configuration of the 6j-symbol. For a general configuration 
one needs to take the square of the symbol. We will not deal with this case. 

We now consider pseudo-isoceles 6j-symbols, which admit an integral formulation. It makes it easy to derive a new 
recurrence relation. We study its relation to the Biedenharn-EUiott identity, and emphasize the underlying geometric 
content. 



A. Recurrence relations for the Isoceles {6j}-symbol 

The so-called "isoceles" 6j-symbols, relevant to the calculation of graviton-like correlation in 3d quantum gravity 
[H, H^j admits a simple integral formulation, 

F{j,k,,lK) EE J dgdhx,{9)Xk{h)xj{9h)xK{gh-') = (-1)'^' '^j k] ^ 

where Xjig) = smdjO / smO is the character of the SU(2) representation of spin j evaluated at the group element g 
with class angle 9. Here dj = (2j -I- 1) is the dimension of the representation. Notice that (— )'^-' = (— )^'° because of 
the parity condition on the two triplets of representations {j + J + K) G N and {k + J + K) e N. 

Thanks to the symmetries of the integrand, the two integrals over SU(2) can be reduced to a single triple integral. 
This is most conveniently written in terms of four angles subjected to a constraint. Let us call a and /3 the class 
angles of g and /i, and (jy^ those of the products gh^. Then we can write ^ as an integral over these four angles with 
the constraint 

cos a cos /3 = -(cos^"'' -I- cos(/)~). (10) 
Observe that pU]) can be written in terms of the characters of the fundamental representation as 

C{g, h) = XL {g)xi (h) - XI (gh) + XI {gh-') = 0. (ii) 

Inserting C{g, h) in the integral ([9]) gives automatically zero. On the other hand, using the recoupling formula 

Xi(5)Xi(5) = Xj_i(.g) + Xj+i(5) (12) 

we can split the integral into a sum of different isosceles 6j-symbols. This method allows us to immediately derive the 
following recurrence relation. 



^ The model can be made less trivial considering the Turaev-Viro model, because a non-vanishing cosmological constant introduces an 
extra quadratic dependence in the Wheeler-DcWitt equation. However we will not pursue this line further. 
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We would like to point out the simplicity with which we derived this relation. As mentioned above, all properties 
of the 6j-symbol are encoded in the Biedenharn-EUiott identity, thus also this one can be derived from it, as we show 
below. However it is remarkable to notice that having an integral representation at disposal paves the way to natural 
and straighforward derivations of recurrence relations. 

Unlike the standard recurrence relations found in the literature, this one has trivial coefficients and all the spins 
are shifted. These features arc achieved at the price of increasing the number of the terms in the relation, from the 
usual three or four to eight. As a difference equation it is second order, since we are dealing with the sum and not 
the difference of terms involving variations of the same variable. Finally, notice that with the help of an additional 
iteration it can be put in a form involving spin shifts of 1 instead of 1/2. Such spin shifts by integer step might be 
useful in some contexts, in particular they respect the parity constraints on the representations. 

In the perspective discussed above, it would be interesting to study whether this equation can be compared with 
a discretization of the Wheeler-De Witt equation on the whole tetrahedron, and we leave this issue open for further 
work. Let us nonetheless point out that even if we take all spins equal to j, (|13p does not give rise to a second order 
equation in j. Therefore also taking (jl3p as a starting point, a homogeneous equation like ([8]) can only be obtained 
approximately. 



B. Relation to the Biedenharn-Elliott identity 

As we pointed out above, all properties of the symbol are captured by the Biedenharn-Elliott identity, including 
the standard recurrence relations such as Therefore, also the new (|13|) should be derived from it. This is indeed 
the case, as we now show. The derivation is a bit lengthy, but it endows (|13[) with a simple geometric meaning. 

The Biedenharn-Elliott identity reads: 

where S is the sum of the nine fixed spins. Geometrically, the 6j-symbol |^ ^ ^ ''^^ seen as a tetrahedron whose 

edges are labelled by the representations, with lengths Ij = J + ^, and its triangles are {j, h, g), [j, a, 6), {g, fc, a) and 
(fc, h, b). The property thus says that the amplitude associated with three tetrahcdra sharing the edge I equals 
that of two tetrahedra glued along the triangle (g, h, j): this exactly states the invariance of the Ponzano-Regge model 
under the 2-3 Pachner move. 

To recover the recurrence relation from this identity, we first specialize / = 1/2. As a consequence, d = g + f3 
and c = h + a for a, = ±i, and the sum over I reduces to two terms, I = k ± ^: 

The geometrical version of this relation is illustrated in Fig. [TJ The starting point of the move corresponds to adding 
a flattened tetrahedron to an initial one, so that the overall result can be seen as a "displacement" of the point D to 
D'. 

Consider the first two 6j-symbols on the right hand side: they correspond to the the tetrahedra ABCD' and 
ACDD', glued by the triangle ACD' with edge spins {a, g + (3,1 = k ±1/2). We then apply the 2-3 move again, under 
the form but now on the new tetrahedron ABCD'. Geometrically, a flattened tetrahedron ACC'D' is glued 

along the triangle ACD' , and the 2-3 move transforms ABCD' into ABC'D', inducing shifts for the links meeting at 
C . We choose the shifts of the links AC and CD' to be a' = a and f3' = — /3. The effect of this choice is to eliminate 
the shift of the spin g due to the first move, while a is shifted by a, exactly like the opposite link h. Finally, each of 
the two terms of the r.h.s. of pS]) gives two terms, BC becoming BC with the two possibilities j 

These two successive moves lead to a relation between the five symbols: 

[ill]- 

and for a fixed a. The coefficients of these relations are generically complicated. However, they become much simpler 
for the pseudo-isoceles case considered in the previous section. Take h = a and g ~ b. Also redefine a — » a — a. The 
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FIG. 1: The identity (|15|l . On the left hand side, the tetrahedron ABCD with edges {g,h, j, k, a,b) is glued via the triangle 
{g,h,j) to the flattened tetrahedron BCDD' with an edge spin being 1/2. On the right hand side, the three tetrahedra ACDD' , 
ABDD' and ABCD' share the edge AD' with spin ± i. The tetrahedra ACDD' and ABDD' are "flattened" because of 
the spin i along DD' . Thus, from the point of view of the initial tetrahedron ABCD, the move can be seen as a "small" 
displacement of the point D to D' , resulting in some small length shifts. 



relations become for a ~ 



and: 



k a + ^ b j ^ 2'^ 2'\k-^abj^ 2' 2' \ k - ^ a b 

+ + + 1)0+,.-. I J} + („ + *-. + i)o + .-„+i){.m 1 1] (18) 

By symmetry, we get two similar relations with a shifted b instead of a. Now, the sum of these four relations 
precisely gives the recurrence relation (fT3|) . Notice first that the left hand sides of the above relations form exactly 
half of the relation, while the right hand sides involve the expected other terms. It then turns out that the sum, 
symmetrized in the exchange (j <-> fc), indeed leads to trivial coefficients. 

We have finally the following picture. Remember that the Biedenharn-EUiott identity is the key to the topological 
invariance of the Ponzano-Reggc model. Equivalently, we can say that it generates and asks for the invariance of 
the vertex amplitude under the 2-3 Pachner move, i.e. under both the diffeomorphism and Hamiltonian constraints 
of Riemannian 3d gravity, since the Ponzano-Regge spin foam model implements these constraints on spin network 
states. Equations (jl7[) and (|18p thus generate and ask for the invariance under two successive moves, which deform 
the original tetrahedron by gluing some flattened tetrahedra. The relation (|13p , which is weaker since it results from a 
sum of (jl7p and (|18p . then implements the invariance under a combination of these moves, allowing for compensation 
between them. 



C. Recurrences from Holonomy Operators 

Another way to derive recurrence relations on the 6j-symbols is to "come back" to spin network functionals and 
look at the action of holonomy operators. Let us consider the spin network functional based on the tetrahedral graph 
labeled by representations ji, .., je^ 

6 4 
(=1 n=l 

Here the intertwiners iy are given by the normalized Wigner's 3j'TO-symbols, and its evaluation at the identity gives 
the 6j-symbol, 
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FIG. 2: The functional <(5{j,}(pi--56) can be represented by a tetrahedron whose links are labelled by the spins ji and the group 
elements gi {i = 1, ■ ■ ■ ,6). The action of the holonomy operator Xi(545556) translates, upon recoupling, into the functional 
depicted above, where intersections of links stand for 3jm-symbols, and the spins k\,k2, fca are summed over. The dashed lines 
are the links carrying the spin j. 



The underlying tetrahedral graph is such that the edges carrying j4, 75 and jg form a triangle. Consider 
now the holonomy operator (545556 )j which simply acts by multiplication on the spin network functional, 
Xj(545556) '/'{j",}(5ii •■, .96)- We can recouple the matrix elements of XjidiQ^de) with the matrix elements of 174, 55 
and 36 already present in the functional ip^j^-j{gi, .., .gg)- This gives a recurrence relation between spin network states, 

X,(545556)^te}(5i,---,56) = E 

which can be represented as in picture^ The coefficients 0?^ = (2/s+l) are the dimensions of the spin k representations. 
Evaluating both sides at = I, we find 



Ji J5 J6 



k4 J4 j] (h J5 j] (h je j 1 f ji j2 js 
je h j2 J 1 j4 h ja J [ js k^ ji J 1 fc4 ke 



(22) 



Such recurrence relations for the 6j-symbol are well-known What is interesting about them is the fact that they 
can be derived by the action of holonomy operators, and thus can be directly extended to more general spin network 
states. Furthermore, the technique then gives recurrence relations to more general {3nj}-symbols, as we will discuss 
below in section [Till In fact, we will see that recurrence relations for the {3rij}-symbols are naturally associated to 
cycles of the underlying graph, as one can guess from the above example. 

But before moving to more general symbols, there is an interesting interpretation of (|22p. Recall that the topological 
invariance of the Ponzano-Regge model is proved showing the invariance of the partition function under the two 
tridimensional Pachner moves, the 2-3 and the 1-4. Invariance under the 2-3 move is a direct consequence of the BE 
identity. Invariance under the 1-4 on the other hand is more subtle. Starting from BE and the orthogonality relation, 
one obtains indeed a formula formally matching one 6j symbol to four summed over four common spins: 



j2 \ j ji j2 js 1 / ( VI ( V'^ \ J -^1 -^2 js 1 j Je js ji 1 f j4 j2 k 1 j h jb J4 



I1...I4 

(23) 

which has the right geometric structure to represent the 1-4 move. However this identity is formal because both 
sides diverge. This is usually cured introducing a cut-off on the spins [3l|, and the resulting model is argued to 
be triangulation independent in the limit in which the cut-off is removed (see [1^ for a discussion). Notice then 
that the cut-offed 1-4 move can be reconstructed from (j22p simply multiplying by dj and summing over j up to the 
desired cut-off. As we send the cut-off to infinity, the extra sum in the right hand side diverges, and the left hand 
side reproduces the correct diverging factor 6{T) = dlj. The recurrence relation ([^^ thus provides an alternative 
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regularization of the 1-4 move, where instead of a cut-off on the spins, we fixed to the vahie j one of the spins being 
summed over. This alternative regularization can be seen as a partial gauge-fixing in the spirit of [s^ . and we will use 
a similar method to gauge-fix the Pachner moves in 4d. Conversely, we may see (|22p as the equality of the summands 
of the above formal 1-4 move for any fixed spin I4. In this sense, the recurrence relations coming from the 1-4 move 
are obtained by gauge-fixing the value of the spin responsible for the divergent factor. 

Furthermore, the relation between (j22[) and (j23[) can be extended to spin network functionals. Indeed, instead of 
acting on the state with a single character, let us act with the delta distribution 6{g4gQgQ) = djXj^g^g^ge). 

The result can be easily obtained from (pij) : 

3(949596) n,,}i9i, ... ,56) = E (-l)^'^+^'^+^-^+^'^+^'^+^'^+'^+'-^+'^+^'4,4..4erf, 

fee n]\H h J3/U h j,jnn,n.sMMM}i9i,---,96), (24) 

The evaluation of this equation at the identity reproduces the formal 1-4 move (|23p . in which the presence of the 
divergent factor 6(1) now becomes obvious. In addition, notice that (5(54(?556) imposes the classical flatness constraint 
F = on the triangle which carries 174, 95 and g^. Therefore, the 1-4 move is really generated through the action of 
the flatness constraint on a spin network state, a further connection between the classical symmetries and quantum 
recurrence relations which was already pointed out in [ll[. 



III. RECURRENCE RELATIONS FOR {3nj}-SYMBOLS 

In the previous section we described three different methods to obtain recurrence relations for the 6j-symbols, (i) 
inserting an existing constraint in the integral representation, (ii) using the defining BE identity, and (iii) acting with 
holonomy operators. With the exception of (|13p . the relations we found are well-known in the literature. When moving 
to higher 3rtj-symbols the literature is rather scarce, with the notable exception of the 9j symbol, see [l3|. Because 
higher 3nj-symbols are relevant in models of quantum gravity, it is useful to show how the techniques described earlier 
can be applied to obtain recurrence relations for general symbols. In this and the next Sections we describe how one 
can use the techniques learned above to generate recurrence relations for general symbols. Below in section |V] we will 
show an application of (i) to the lOj symbol. 

Consider an arbitrary {3nj}-symbol, represented by a closed graph made of 3-valent vertices and whose links carry 
SU(2) irreducible representations, collectively denoted {ji}. The evaluation of the graph can be obtained according to 
conventional rules assigning Clebsch-Gordan coefficients (or 3TOj-coefficients) to vertices. To unambiguously represent 
symbols with such spin network graphs, we use the conventions of (l7| for the orientation of links and vertices. 

A straighforward way to obtain recurrence relations is to exploit the decomposition of the 3nj'-symbols into lower 
symbols. This decomposition can be made isolating a part of the graph containing a cycle, as shown in Fig 2. In the 
simplest case, the cycle contains only three links, and the symbol decomposes as {3nj} = {6j}{3(n — l)i}. Foi' a 
cycle with four links, {3nj} = ^i{9j}i{3{n — l)j}i, and so on. 

We show how certain types of recurrence relations can be obtained using these decompositions and applying the 
known relations of the 6j and the {9j}. This procedure works quite well for decompositions built from 3-cycles and 
4-cycles, but becomes cumbersome for larger cycles. Then, in the next Section we focus on 71 = 5, the relevant case to 
4d models of quantum gravity, and describe how these recurrence relations can be obtained from a regularized version 
of the 2-4 move. Interestingly, this derivation will open a window on new recurrence relations, which apply to graphs 
with larger cycles. 



A. Recurrence relations from smaller symbols 



The simplest situation is a {3nj} symbol in which we can isolate a cycle formed by three links, see top panel of 
Figure[3l In this case the symbol is reducible, as it can be written as a smaller {3(n — l)j}-symbol times a 6j-symbol. 
Since the labels of the 6j-symbol do not enter the larger symbol, we can trivially use the recurrence relations of the 
6j to infer relations on the {3nj} symbol, for instance 

h h h] (h h I3] ( h h h' 

A+i{j) {31 + I 32 33 } + Aoij) { 31 32 h ) + ^-i(j) < 31 - 1 32 33 } = 0. (25) 
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FIG. 3: The top figure displays a reducible symbol: the cycle (jij^js) can be factorized because there is an unique intertwiner 
from Zi dg) ?2 <8) h to C. One can then use the recurrence relations of the 6j-symbol to shift the spins (ji). For a cycle made of 
four links, one needs to sum over the intertwiners z between I41S1I2 ^ h One can then use recurrence relations which shift 
some (ji) provided they do not depend on z, but only on the spins which join the cycle. 



Consider next a symbol containing a cycle made of four links, bottom panel of Figure [3l We can still apply the 
recoupling theory to write the symbol in terms of smaller ones, but this time a summation over smaller symbols 
appears: the original symbol is not reducible this way. This docs not stop us from applying recurrence relations to 
the smaller symbols, but care is needed to avoid including into them the spin being summed over. Using the example 
in the figure, one can derive relations for the symbol to the left starting from relations for the {9j}-symbol on the 
right which do not include the spin z intertwining the representations {h} together. For instance. 



in + I2+ h + 1)02 + j3 - l2){-jA + J3 + /3)04 + ^3 ^ J3 + 1) 



Jl + \ k i4 + I ' 
2 '2 J3 - 2 



.72 



+ 



ijl+h+j2 + 2){ji-h+j2 + l) 
(j2 - J3 + h){-]2 + J3 + h + l){k + J3 + ^3 + 2)(j4 + J3 - ^3 + 1) 



(jl+?l+j2 + 2)(ji-;i+j2 + l) 
{h +I2+ .13 + 2) {j2 + J3 -I2 + 1) (j4 + i3 + ^3 + 2) (j4 + ^3 - ^3 + 1) 



(jl+il-i2 + l)(-il+^l+j2) 

(.72 - k +h + 1)(-J2 + k + ^2)(-i4 + ]3 + ^3)^4 + ^3 - 33 + 1) 
(jl +^1 -i2 + l)(-il +^1 + j2) 

(.n+'4+j4+2)0-i+/4-j4 + l) 



jl + \ li ji + \ 

j2-\ h k + \ 

' h + \ h .14 + ^'^ 

h z h 
h + ^ h k + \ 



= (2j2 + l)(2i3 + l) 



- [jl +h- 32 + 1) (-71 +h+ 32 ) (.71 + ^1 + .72 + 2) {ji -li+32 + l) 
Other relations can be obtained from other relations of the {9j}-symbol given in [iTj . 
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Notice that the coefficients of recurrence relations obtained this way will depend on the spins of the cycle and on 
those of the links directly joining the cycle. This property also holds for cycles made of five and six links. However the 
procedure can become quite cumbersome for such cycles. Below we will derive an alternative which is more fruitful 
for large cycles. 

Let us also mention an alternative used by Jang p^ . A {3(n + l)j {-symbol can be expressed in terms of a sum 
of products of a {3rij}-symbol and typically some 6j-symbols. Then, using symmetries of the {3(n + l)j}-symbol, 
one can get non-trivial relations between the {Bnjj-symbols, with some 6j as coefficients. However, it seems that this 
method works for symbols of the first kind but may not be efficient for other kinds of symbols. Moreover, it does not 
provide any geometric picture such as simplicial moves, in contrast to the method we will use. 

IV. THE 15J SYMBOL AND 4D BF THEORY 

Let us now fix n = 5. The resulting 15j symbol is not unique, unlike the 6j and 9j symbols. Using the terminology 
of [13, [3 J one can distinguish five different irreducible kinds. Depending on the kind, different types of cycles will be 
present, and one can derive relations for the 15j symbol proceeding as described above. All links of the 15j-symbols 
of the first and second kinds can be included in cycles made of four links, but five links are necessary for the 15j of 
the fifth kind. The symbols of the third and fourth kinds contain both types of cycles. 

The 15j symbol was used by Ooguri [T^ as the vertex amplitude of four-dimensional BF theory with gauge group 
SU(2), a theory often used as a starting point for spin foam models of quantum gravity. To see how the 15 symbol 
enters the game, let us recall that Ooguri's model is built triangulating spacetime with a collection of flat 4-simplices 
glued together, and assigning an SU(2) element gt to each tetrahedron t. Within the triangulation, each triangle / 
will be shared by many tetrahedra, and one can define the quantity gj = OtD/fft; called the holonomy associated to 
each /. Finally, the partition function of the model is the integral over all the gt of Y[f S{gf)- These delta functions 
can be seen as imposing a discretized version of the flatness equation which defines the topological BF field theory. 

The integrals can be evaluated using the Planchercl formula 6{gf) = X^jg^ d,jXj{gf) at each triangle and then 
writing Xjigf) = XjiY[fDf9t) ^^e trace of a product of matrices, one for each tetrahedron sharing the triangle. 
Noticing that each tetrahedron contains four triangles, there is only four matrices containing a given gt in the partition 
function, which then factorizes into the following integrals, 

r ^ 

/ dgt ^D^if\gt)^ idi,,^. (27) 

JSU(2) 

Here idinv is the identity on the invariant subspace of the tensor product of the four representations jf meeting 
at t. An orthonormal basis of such 4-valent intertwiners is given by a tree expansion, i.e. an expansion onto 3- 
valent intertwiners: one has to choose one of the three possible pairings between the four jf, and sum over the 
internal representation. For instance, an intertwiner is |(jij2)j (j3j4); where i satisfies: |ji — J2I < « < ji + j2 and 
|j3 - ji\ < i < ja + J4- The identity is given by: 

idinv = X! {h3i);i) ((jij2), {jzjAYA (28) 

i 

The contractions of such intertwiners at each 4-simplcx give rise to a 15j-symbol. With respect to the original 
triangulation, the model naturally assigns a spin jf to each triangle, an intertwiner it to each tetrahedron, and the 
resulting 15j symbol as the amplitude of a 4-simplcx. 

The attentive reader will have noticed that which 15-symbol emerges depends on the choice of pairing in (|28p for 
each of the five intertwiners in the 4-simplex. It is possible to write the model in terms of a single kind of 15j-symbol, 
but generally this requires using different pairings for each of the two intertwiners entering (|28p . Gauge invariance 
is guaranteed by the subsequent emergence of appropriate 6j-symbols at the tetrahedra recoupling the two different 
pairings.^ 

A natural question is whether exploring all the possible pairings one can reconstruct all the five different kinds of 
15 symbols. It is easy to show that one can indeed recover the 15j of the first, third, fourth and fifth kinds (in the 
terminology of [isj). along with reducible symbols. On the other hand, we were not able to prove that the second 



^ This is often overviewed in the Uterature on spin foam models, although these 6j-symbols are part of the model originally defined by 
Ooguri. 
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FIG. 4; The left picture describes the configuration with two 4-simplices, while the right picture displays the structure of four 
4-simplices which are glued to each other along different tetrahedra. The arrows give the orientations of the dual edges and 
the boxes stand for integration of the group elements carried by these dual edges. Once the right hand side is regularised, 
the two situations have the same BF amplitude. To regularise it, we have dropped the flatness condition for the face which 
passes through I, q and p. As a consequence, one obtains only one 15j-symbol on the right hand side, since only three lines 
pass through I, q and p: three 4-simplex amplitudes are reduced to {12j}-symbols. 

kind can also be obtained in this way. We do not know whether 15j-symbols of the second kind can appear in the 
spin foam model. 

The interest in Ooguri's model model is its topological invariance. which gene ralizes to four dimensions the analogue 
property of the Ponzano-Regge model. The proof uses again Pachner moves U^^. In4d, there are three independent 
moves, which are discussed in details in [30|: 3-3, 2-4 and 1-5. When applied to the vertex amplitude of the model, 
the 15 symbol, both the second and the third moves are divergent [s^l. As for the 1-4 move of the 6j, the source of 
the divergences arises from redundant delta functions imposing SU(2) flatness in the bulk. 

Let us now come back to the recurrence relations that we can construct using the strategy of the previous section. 
Since the 15j symbol gives rise to a topological model, one might wonder whether these relations are again a conse- 
quence of a more fundamental identity, like the BE identity for the 6j, and whether like the BE identity, they are 
related to a Pachner move. We now show that this is indeed the case: the recurrence relations are a consequence of 
the (regularized) 2-4 move. Since the 2-4 move is a component of the invariance under homeomorphisms of the model, 
the recurrence relations can be seen as difference equations contributing to the implementation of this symmetry, i.e. 
as part of a discrete version of the classical constraints of BF theory. 

However the relevant move is divergent, so we need to regularize it first. We do so in the the subsection, and then 
move to the derivation of recurrence relations from it. 



A. Regularizing the 2-4 move 

Since the model is intrinsically defined without choices of pairings, but rather in terms of group integrations, it 
is convenient to consider the move in the group picture, see figure ID Black boxes stand both for edges dual to 
tetrahedra and integrations over the corresponding group elements gt, while lines denote the boundaries of the faces 
dual to triangles. Each box is crossed by four lines. On the two sides of the move, there are eight boundary tetrahedra, 
with elements ((7i)i=i,....8- On the left side, the two 4-simplices are glued along a single tetrahedron, with element g. 
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Using the orientations of dual edges given in figure [H the amphtude of the left hand side is 

8 

/ H dgt D^=^-\g^g^^) D^^^-\g^g^^) D^^^-\g,g^') D^^^^^g^g^') D^^^^\g,g^') D^^^--\g,g^') D^^^^^g^g^') 
t=i 

D^'^^Kg^ge') D^'''\g^g^') D^'"'\g^g^') D^'^^^g^g^') D^^-^\g,g,') D^^^^Hg^g,') 

J dg (.91 335"') D^'''^ (.92 ,9 Se"') D^'''\g3 5 57"') D^'''\g4 g g^'), (29) 

where tensor products are intended between the matrices Z?^^-* . After integration over the elements gt for the eight 
boundary tetrahedra, one can contract the free indices of intertwiners, corresponding to the free ends of the open 
lines, with eight chosen intertwiners. This process includes choosing pairings of the virtual spins. The integration 
over g, i.e. over the tetrahedron shared by the two 4-simplices, produces the sum (|28|) . and this leads precisely to the 
left picture of figure [S) 

On the right hand side, one has to integrate over six elements corresponding to the six bulk tetrahedra. Open 
lines correspond to boundary triangles, closed lines correspond to triangles in the bulk and are shared by exactly 
three tetrahedra (and three 4-simplices). Each bulk triangle contributes a factor S{Y[f-^j gt) to the amplitude. By 
inspection, one of this delta is redundant, and this causes the (5(1) divergence, exactly like in the 1-4 move of 3d 
gravity. This can be removed brutally dropping one of the four flatness conditions, say S{lqp~^), since the constraint 
Iqp^^ = I is ensured by the other delta functions. The resulting finite amplitude is 

/ ]ldgtD^'^^\gig^')D^'^-Hg2g,')D^'^^\g3gj')D^^-^\g,g,') f dgdhdkdldpdq 
•' t=i 

S{glh-')S{gpk~')5{hqk-') D^^^-\g, g g^') D<-^^^<^\g, g g^') D<^^^^\g, h g^') D^^^^\g, h g^') 
D^^--\g,kg^^) D^^^^\g,kg^^) D^'-^Hg^l g^') D^^^^^Hg^l g^') 

D^'''\g2Pgl^) D'-^<^'\gepgs') D^^^'^g^qg^') D^^^'^gr q g^'). (30) 

The bulk integrals can be performed as usual, expanding the delta functions with the Plancherel formula and using 
repeatedly (P7|) . Once again, the latter step introduces a choice of pairing. Notice however that the result will involve 
only one 15j symbol: because a delta function has been removed on one triangle, the group elements Z, p and q now 
appear only three times. Thus the amplitude for the three 4-simpliccs sharing this triangle arc reduced to smaller 
12j-symbols. 

The two expressions coincide after the regularization, since all we have changed was to remove a redundant delta 
function. It can also be explicitly proved using the graphical method of The scheme is basically as follows. 

Notice first that g, on the left hand side, can be reabsorbed on the right of 31, 32, 33 and 34, using the translation 
invariancc of the Haar measures. We briefly describe how the right hand side reduces to the same amplitude. The 
same method can be used to eliminate h and q for instance, so that k is forced to be the identity, g can also be 
reabsorbed, enforcing I = id. Finally, integrating p is trivial because one gets from the previous simpliflcations S{p). 

Because the divergence is exactly the same as in the 3d 1-4 move, we can obtain formulas similar to (|22p . Indeed, 
instead of removing the delta function S{lqp~^), let us keep a component, say of spin J, of its Fourier expansion, 
Xj{Wp~^) in the integrand. Performing the integrals and choosing the pairings, the right hand side of the move is 
now made of four 15js, with a dependence on J. Let us see how it changes the amplitude of the left hand side. Since 
the other delta functions (in the right hand side) enforce lqp~^ = I, the character XjQ-1P~^) can be reduced to dj 
which comes in factor of the above expression. Thus, the two 15js of the left hand side are only changed by a factor 
dj, like in 

We now show how this regularized move can be used to obtain the previous recurrences and new ones. 

B. Recurrence relations from the 2-4 move 

Following the derivation of p3|l from the BE identity, we look at the regularized 2-4 move with some of the boundary 
spins fixed. The choice has to be careful, to ensure that we still have one 15j symbol on the right hand side. The 
situation is illustrated by Figure^ we consider a 4-simplex, with points a, 6, c, d, e, and we glue another one, denoted 
(aa'cde), along the tetrahedron (acde). This will be the left hand side of the 2-4 move. Every triangle is colored with 
a fixed SU(2) representation, and every tetrahedron, except (acde), by an fixed intertwiner: these are the 'external', 
or boundary, data. The tetrahedron (acde) being shared by the two 4-simplices, the corresponding intertwiner is 
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FIG. 5: The two sides of the 2-4 move. On the left hand side, the 4-simpUces (abcde) and (aa'cde) are glued along the 
tetrahedron (acde). On the right hand side the four 4-simplices (a'bcde), (a'bacd), {a' bade) and (a'bace). This configuration 
has four 'bulk' triangles (given by (a'bc), (a'bd), (a'be) and (aa'6)), which are shared by three 4-simplices, and six 'bulk' 
tetrahedra (made of the points a and b together with any choice of two others), which are shared by two 4-simplices. When 
the spins of the triangles (aa'e) and (aa'c) are taken to be zero, we may consider a' as being very close to a. Then, the areas of 
the initial 4-simplex (abode) on the left and those of (a'bcde) on the right only difi'er by some slight shifts for the three triangles 
sharing the edge (ad). 

summed over to glue them (due to the integral over g in (|29p ). The right hand side configuration is made of four 
4-simplices respectively obtained by dropping the points a, c, d and e. 

To write the indices of spins and intertwiners, we will use the standard dual notation: wc denote a simplex label by 
only writing the points which are not vertices of the simplex. Typically, the representation coloring the triangle made 
of the points (cde) is jaa'b and that coloring the tetrahedron (bcde) is iaa'- To regularize the move, we need to drop 
the delta function of one internal face of the right hand side: we remove that of the triangle (aa'b). Thus, only the 
amplitude of the 4-simplcx which do not share this triangle, i.e. (a'bcde), is actually a 15], while those of the other 
three 4-simplices are reduced to some smaller 12j-symbols, as indicated in Figure [S] 

We now fix two boundary spins to zero: jbcd — jbde ~ 0. The corresponding links in the spin network graphs are 
represented by dotted and dashed lines in Figure [51 This enforces the equality jatd = ja'bd between boundary spins. 
As a consequence, the 15j of (aa'cde) on the left hand side turns into a (possibly reducible) 9j-symbol. It turns out 
that for any pairing of the boundary intertwiner iat, we can always expand the summed intertwiner ia'b so that this 
9j is further reduced to a 6j-symbol or at least a product of two 6j-symbols. This way, the left hand side contains 
only one 15j and one or two 6js. 

A similar reduction applies to the right hand side. In particular, among the spins which arc summed over, two 
become fixed: jacd ~ ja'cd and jade — ja'de- Also, the three 12j-symbols arc reduced to products of 6j-symbols, whose 
precise forms depend on the choice of pairings for the boundary intertwiners. 

We are then in a position to write relations with one 15j on each side, and with coefficients being 6j-symbols. It 
is also easy to see that all these 6j-symbols depend on the spin A = jbce- This spin also parametrizes the ranges of 
the data which are summed: typically, one gets from Figure [H] that \ja'ce — A| < jace < ja'ce + A. Thus, exactly as 
we specialize one spin to 1/2 in the BE identity to get recurrence relations (or to 1, to get recurrences with shifts 
of 1), wc can specialize A to small values and this way obtain recurrence relations. Notice also that the conditions 
jbcd — jbde = cntail some relations between the boundary spins, which also depend on A, namely: jabc = ja'bc + 
and jabe ~ ja'be + P , whcrc a and (3 are such that: —A < a, /3 < A. 

Before we give the general formulas, we would like to have a geometric picture of what we have done with the 2-4 
move. In 3d, the geometric content of the move and of the recurrence relations in terms of gluing some flattened 
tetrahedra is supported by the interpretation of spins as quantum lengths. In 4d, a natural interpretation of the 
spins labelling triangles and of the intertwiners labelling tetrahedra comes from the geometric quantization of a single 
tetrahedron and the similarity of BF theory and loop quantum gravity at the kinematical level. It invites us to think 
of the spins j as being related to the triangle areas, and the intertwiners i as areas of parallelograms chopping each 
tetrahedron in half (notice that the three possible pairings of each intertwiner correspond to the different ways of 
chopping each tetrahedron in half) (l9| . 

Then, choosing jbcd = jbde = means that the triangles (aa'e) and (aa'c) are "very small", so that the 4-simplex 
(aa'cde) is flattened. From the point of view of the 4-simplex (abcde), the overall result can be seen as a small 
"displacement" of the point a to a' , see Figure[5] On the one hand, the relations jacd = ja'cd and jade = ja'de mean 
that the new triangles (a'be) and (a'bc) have the same areas as the initial triangles (abe) and (abc). On the other hand, 
the areas of the triangles sharing the edge (ad) are shifted when moving the point a: as wc have seen, jabc = ja'bc + ct 
and jabe = ja'be + P, and ja'ce IS also shifted by a quantity depending on A. Varying A corresponds to moving the 
point a' so that the area of the triangle (aa'd) changes, while keeping those of (aa'c) and (aa'e) fixed and very small. 
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FIG. 6: The regularized 4-2 Pachner move. Only one graph on the right hand side is a 15j-symbol since the flatness condition 
for the triangle (aa'b) has been dropped, x stands for the summed colors, whose corresponding links and vertices (intertwiners 
are not displayed) are marked. To get formulas involving only one 15j-symbol on each side of the move, we take jbcd = jtde ~ 0, 
which are represented by dotted and dashed lines, so that many graphs reduce to (products of) 6j-symbols. Looking at small 
values of A = ji,ce enables to get interesting explicit recurrences. 



Whether or not the intertwiners are shifted when moving a to a' (for instance, whether or not iab = ia'b) depends 
on the chosen pairings. More precisely, we look at the tree expansion of an intertwiner in the spin network picture. 
If two shifted spins meet at a common vertex of the intertwiner, then the spin i of the latter will not be changed. 
But if the line carrying the spin i is inserted between two shifted spins, it will also be shifted. This way, one can 
obtain recurrence relations for 15j-symbols whose cycles have different numbers of links, by changing the pairings of 
intertwiners. 

Relations for cycles made of three and four links turn out to be precisely those for 6j and {9j}-symbols, as expected. 
Let us present new relations, adapted to cycles of five and six links. There are two ways of inserting intertwiners 
between the shifted links to get a cycle of five links. Indeed, one can insert, or not, ia'b between ja'bc and ja'be on the 
left hand side, while if ia'b is inserted, it is then equivalent to insert ia'e or ia'c- We thus get two different relations: 



ioti + A ^ . • ~1 r \ ■ ■ ~1 I *a'e Ja'bc 

\R ^ ^ J A *a'fc «afc I J A Ia'b iab 

ia'b = Hab-^\ 



^ ' ^ '-'Aja'bd ]a'bc + a ]a'bc\\jaa'b ia'be+P ia'be^ ^ ^"'"^ 



E7 7 J ^a'c Ja'bc Ja'ce I I ^a'e ^ae 1 J ^a'e ^ae I J - , ri 



jace ja'bc H~ J \ja'de jabe ^~ P ja'be J \jaa'e jace ja 



^ab ^ae Ja'bc H" ^ 



ace 



where a and /3 are half-integers between —A and A, dj = 2j-hl and R — ja'de'^ja'bd~jaa'e~jaa'b — iac + '^ja'bc-^Oi-\-X; 

^ I Ja'bc ^a'e ^a'c i ( \ 

iab Ja'bc Ja'be \\ ■ ■ I {-if d,^J,_d,\ .''^'^'^ 

A Ja'be +P Ja'be + a] | f ^ ^ ^ ^ »■= {ja'de Ja'be +P Ja'bc 



A l-a'p. ^ae 1 I A ^a'r, ^ar, 1 f A 



Ja'bc ~^ ^ ^ae ^ 



ac 



X i/' U . [><{." ^."''^ \ { Ja'be + 13 Jace ) (32) 

\Jaa'e Jace Ja'ce] \Ja'cd Ja'bc \ ^ Ja'bc] \Jaa'c Jace Ja'i 

for S ^iab^ iac - ia'c + '^^a'e + 2ja'bc + « + jade + jaa'e + jacd + jaa'c- 




FIG. 7: This picture represents the 15j-symbols respectively appearing in pif) . (|32|) . which give recurrences for a cycle of five 
links, and in (|33p with a cycle of six links. Notice that only the orientations of the vertices and links of the cycles and of the 
links joining them are required. 

There is only one way to make a cycle of six links, which gives the following relation: 



^ 1 Ja'bd Ja'bc ^ Ja'bc J 1 Jaa'b Ja'be 



'be 



Ja'be Ja'ce Ja'bc 



E 



ja'de ja'be ~^ P ja'be] \jaa'e jaee ja't 



-1) di^^dj^^^di_ 



i ^ ^a'c ^ac 11 ^ ^a'c ^ac 11 ■ i /3 ■ ■ i I /ooA 

X i . , ^ , , , , M , , , , , M 3a'be + /? Jace Ja'bc + « } (33) 

Ja'ea Ja'bc ^ ^ Ja'bc I [Jaa'e Jace Ja'ce 



^ab ^ae ^a 



with T 2ja'hc + 2\ + a + ia'e + «ac + ja'bd + .7aa'6 + jaa'c + jaa'e + ja'de + ja'cd- The cocfScicnts of this cquation are 
naturally those of the left hand side of (pij) with those of the right hand side of ((5^ . The 15j-synibols described by 
([?T|) , ([5^ and ((551) represented in figure |7| 

These relations involve 15j-symbols whose spins are shifted by half-integers between —A and A. To get explicit 
recurrences, the simplest choice is to take A = i. Then, the free parameters are a,(3 = ±i, and every spin of the 
cycles is shifted. In each of three above formulas, exactly three triangle spins are changed, since the cycles only differ 
from the insertion of intertwining spins (which are also changed). The choice A = 1 is also of interest. In this case, 
one is free to choose a, (3 — ±1. The case a, (3 — Q simplifies the equations since it enables to leave invariant two 
spins for each cycle. Geometrically, it means that only one triangle area, ja'ce, is now changed by the move, together 
with two or three intertwiners. These special cases are summed up in Figure HI 

In the above formulas, only the spins along a cycle of the spin network graphs are shifted. We emphasize that the 
coefScients of these relations only depend on the spins of the cycles, and on those of the links joining them. As a 
consequence, they apply to any invariant {3ri,j}-symbol containing such cycles. This gives us a generalization of the 
method described above for smaller symbols to symbols containing cycles of 5 and 6 links. 

Summarizing, we started from the regularized 2-4 move, and obtained recurrence relations by choosing some specific 
boundary spins, exactly like in 3d with the BE identity. We could also have followed another programme proposed in 
3d: by acting with holonomy operators on spin network functional. We used this method to provide a well-defined 
expression of the 1-4 move and at the same time get some recurrence relations. The same can be done in 4d to provide 
a definition of the regularized 1-5 move. In this case, the naive move diverges due to a factor <5(I)^, which means that 
four delta functions are redundant. Choosing a specific spin ji among the Fourier expansion of each of these delta 
only contributes to the left hand side by a multiplicative factor, like for (|22l) . Moreover it amounts to acting on the 
spin network functional associated to the 15j-symbol with the characters in the spin ji representations, Y\a=i Xjii^i) 
if the elements Gi are the holonomies around the four divergent faces. As we discussed in Section Hi CI the divergent 
move can be reconstructed by acting with Y[i=i ^i^i)-! leading to a precise relation between the quantization of the 
classical constraint F = and the 1-5 move. 

We can also get this way non-trivial relations between 15j-symbols by simply fixing the four spins ji. However, 
invariance of Ooguri's model under the 1-5 move follows from that under the 2-4 move (with orthogonalities relations) 
- like in 3d, the BE identity together with an orthogonality relation can be used to prove the 1-4 invariance. Thus, the 
recurrence relations we would get from the 1-5 move can also be derived from the above and more general formulas. 
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FIG. 8: This array sums up the resuhs obtained for cycles made of five and six links, with arguments changed by | or 1. For 
simplicity, the coefficients of the relations, consisting in 6j-symbols which can be explicitly evaluated, have been omitted. The 
recurrence relations act on any 5-cycles or 6-cycles of any SU(2) invariant graph. For A = i, one can choose the values of the 
free parameters a, (3 — ±i, and all the links of the cycle are shifted. For A = 1, the choice a — f3 = enables to simplify the 
structure of the equation by leaving some spins invariant. 



Although our recurrence formulas can be applied to any invariant symbol, the geometric interpretation coming from 
the 2-4 move only holds for 15j-symbols, due to their relation to 4d topological BF theory. Nevertheless, we may 
expect part of this interpretation to also hold for spin foam models aiming at describing 4d quantum gravity. Indeed, 
quantum states in loop quantum gravity are built with SU(2) spin network functionals, whose spins and intertwiners 
are eigenvalues of area and volume operators. If a spin foam model provides the dynamics of these spin network 
states, we may then expect that recurrence relations for the spin foam amplitudes translate the classical constraints 
(i.e. symmetries) of the theory at the quantum level. Thus, the shifted spins would correspond to shifts of areas, and 
it is reasonable to think of describing such shifts in terms of elementary geometric moves such as the one presented 
here. 



V. A RECURRENCE RELATION FOR THE BC VERTEX 



In this section we consider another special 3rij-symbol singled out by spin foam models of quantum gravity, the 
lOj symbol used in the BC vertex amplitude. This admits an integral representation like the isosceles 6j, and we can 
derive a recurrence relation for it using method (i) of inserting a constraint in the integral. This constraint is related 
to a basic geometric property, namely the closure of the 4-simplcx. We also propose an interpretation of the different 
terms entering the relation using the moves studied above. 
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A. From the Gram matrix constraint to the recurrence relation 



We start with the expression of the Barrett-Crane lOj-symbol for a Euchdcan 4-simplex as an integral over five 
copies of the SU(2) group as given in [H [13, [H, [Ifi] : 

{m = j \{dga X{x,A9a^9b)- (34) 

a=l a<b 

It is a function of the ten representations jab G N/2 labehng the ten triangles of the 4-simplex. The characters Xj{9) 
are defined as the trace of the group element in the j-representation of SU(2) and depends entirely on its class angle 
(half the rotation angle): 

XA9) = ^ = U,,{cos6). (35) 

The dimension of the ^-representation is given by dj = (2j -I- 1) S N. The ?7„ are the Chebyshcv polynomials of the 
second kind and can be defined through the following generating functional: 

E^«(-)^"-T^^- (36) 

riGN 

The most useful expression of the lOj-symbol for both analytical and numerical computations is in term of the ten 
class angles of the group elements ga^9b- Since these group elements are obviously not independent, this results in a 
non-trivial measure on the ten class angles dab- In fact, the measure is given simply by a constraint [25| : 

{lOj} = 4 / [d0r6{dctG[e]) n sin(2jab + l)0ab. (37) 

a<b 

The Gram matrix is a symmetric 5x5 matrix defined as Gab = cos 9ab with the obvious rule 9aa = 0, cos Baa = 1- This 
constraint contains all the geometric information and allows to relate the lOj-symbol to the Regge amplitude 

for a geometric Riemannian 4-simplex in the large j asymptotics (23 . [25l . [26j . We can rewrite the previous integral 
formula in term of the Chebyshev polynomials by defining the ten variables Xab = cos 9ab '■ 

im = 4 r\dxrd{g[x])l[U2Uxab), (38) 

•^-1 a<b 

with the new symmetric Gram matrix Qab = Xab for off-diagonal elements and Qaa = 1 on the diagonal. 

We would like to use the recurrence relation on the Chebyshev polynomials to build a recurrence relation for 
lOj-symbol: 

U2] + l{x) = 2xU2j{x) - U2J-1- 

This relation comes from the simple trigonometric formula: 

sin(2j + 2)6 + sm{2j)9 = 2 cos 6* sin(2j + 1)6*. (39) 

The idea is to apply this formula to the integral expression for the lOj-symbol ((37|) and to combine the cos 9 factors 
so to form the determinant of the Gram matrix det(G). Indeed the determinant reads as a sum over all permutations 
between five elements: 



det(G) = ^e{a)Y[cos9aaia), 

(7 a—1 

where e(a) is the signature of the permutation. We write the following recurrence relation: 

= / 5{Act G[9])Act{G) \{ sin(2j„fc + l)9ab. (40) 



a<b 



For each cos 9 factor coming the determinant, we get a linear combination of an upper and lower shift on the corre- 
sponding representation label using equation (j39p . 
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The trivial permutation a ~ I has no cos^ factors since Va,I(a) = a; thus it gives back the original lOj-symbol 
with no shift on the representation labels. There are no permutations fixing only four elements. The next step are 
permutations fixing only three elements, say a = 1,2,3. Since it permutes the elements a = 4, 5, we get a factor 
— cos^ 64^,: 

- cos^ sin(2j + 1)61 = -i sin(2j + 1)6* - i (sin(2j + 3)9 + sin(2j - 1)61) . 

We get a linear combination of the original lOj-symbol together with other lOj-symbols where the representation label 
j45 is shifted by ±1. Then come the permutations fixing only two elements, say a = 1, 2. The action on the remaining 
elements a = 3, 4, 5 can be either of the two cyclic permutations on three elements. These two permutations give 
the same factor + cos ^34 cos ^45 cos ^35 . This leads to a linear combination of upper and lower shifts by ±i on the 
representation labels ^34, j45 and jas. For permutations fixing only one element, say a = 1, we get either have a cyclic 
permutations of the four remaining elements (6 possibilities) or two permutations on the two pairs (3 possibilities). 
The 4-cycles will give factors of the type + cos 6*23 cos 6*34 cos 6^45 cos 6*25 which will leads to ±i shifts on the relevant 
representation labels. On the other hand, the two pairs of permutations give factors of the type cos^ 6*23 cos^ 6'45 which 
leads to ±1 shifts on the corresponding representation labels. Finally, we have to look at permutations that do not 
leave any of the five elements invariant. Either, we have a 5-cycle (4! possibilities) or a combination of a 3-cycle and 
a 2-cycle (10 x 2 possibihties). At the end of the day, we check that: 1 + 10 + (2 x 10) + 5 x (6 + 3) + (4! + 20) = 5!. 

This procedure leads to a recurrence relation on the lOj-symbols relating the original amplitude with spins jab, 
amplitudes with up to five spins shifted by ±i and/or up to two spins shifted by ±1 (due to 2-cycles). The complete 
expression involves 5! — 120 terms, so we do not write it here. On the one hand, such a recurrence relation could be 
used to investigate the asymptotics and the possible topological properties of the lOj-symbol or more generally the 
behavior of (graviton) correlations in the Barrett-Crane(-like) models. However, on the other hand, we are aware that 
the resulting formula can not yet be used to numerically compute the lOj-symbols by recurrence since it does not yet 
express a single lOj-symbol in terms of symbols with lower spins. Under the present form, the practical use of this 
recurrence relation for calculation purposes would require an infinite number of initial conditions. But wc do hope to 
find some trick to turn it into a recurrence relation useful for actual computations. 

Notice that the recurrence relation obtained in this way is a priori a fourth order difference equation. We will 
comment on this below. 



B. The recurrence as a closure constraint 

The recurrence relation obtained for the BC vertex is quite complicated since it results from the expansion of a 5 x 5 
determinant. So it would be nice to identify a simple quantity satisfying this equation, at least in the asymptotical 
regime. Obvious solutions are given by constants, since: e{a) = 0. We now look for more geometric solutions. 

To this aim, we emphasize that the recurrence relation has been obtained for the lOj-symbol by inserting into its 
integral representation the quantity: det (cos 6*06). Its vanishing is precisely the necessary condition so that the angles 
6ab can be the dihedral angles of a genuine, i.e. closed, 4-simplex. Furthermore, from the ten spins, interpreted as 
the ten areas of a 4-simplex via the formula A f = dj, , one can compute the associated dihedral angles Oab{j) which 
satisfy this condition. In the asymptotics [H, [2^, [2^, the oscillatory part of the lOj-symbol then involves the Regge 
action for the 4-simplex, evaluated in terms of the areas dj^ . Thus a natural candidate to satisfy the recurrence is the 
Regge action for a 4-simplex. More precisely, we will show that the cosine of the Regge action, cos{S{jab) + a), as a 
function of the triangle areas {2 jab + 1): 

S = Y.^2jab + l)eab, (41) 
a<b 

where a an arbitrary phase, satisfies the recurrence in the asymptotics. Furthermore, the key identity to this result 
is indeed the closure of the 4-simplex, encoded into: 

det (cos 6ab) = 0. (42) 

The recurrence relation shifts several spins together by ±i according to the previous discussion. We first imagine 
the case of a shift of only one spin, say ji2- According to the trigonometric identity ((39)) . each shift of \ is accompanied 
by an equivalent term with a shift of — ^ . Also following the analysis of Schulten and Gordon about the asymptotics 
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of the 6j [|, we consider that the spins are large and expand Sr to first order, neglecting its second derivatives. It 
gives: 

i cos(5fl(ji2 + + a) + i cos(5fl,(ji2 =^Y1 (^osi^SnUu) + a + 77^i2(ji2)) (43) 

ri=±l 

cos(^5ij(ji2) +a) cos^iaOah) (44) 

To obtain this result, we have used the Schlafli identity, X]a<h(^Ja6 + ^)Sdab = 0, so that the variation of Sr with 
respect to ji2 is simply 20i2. 

The above expression can be easily extended to the case of an arbitrary number of shifted spins. Let us denote the 
shifted spins jij for some i < j. Then, we obtain the cosine of the Rcgge action times the product of the cosines of 
the angles Oij: 

i cos(^Sii{j,j +T],j^) + = ^ cos(^Sii{jab) + a + Yvzj(^ij{.iab)^ (45) 

Vij=±l »?,j=±l 

= cos (^5*77 (jab) + ]^COS%(jafc) (46) 

Remember now that the shifts in the recurrence we are interested in are not arbitrary but correspond to permutations 
of the links. Moreover, the sign of each term is determined by the signature of the corresponding permutation. We 
thus obtain: 

5 

^e(cr) cos{SR + a) J| cosdaaia) = cos{Sr + a) dctG[e] = 0. (47) 

a a—1 

This obviously holds because the angles 9ab are the dihedral angles determined by the areas. 

This result shows that the recurrence found for the lOj-symbol is related to a simple geometric property of a 
single 4-simplex, namely its closure, unlike the usual recurrence relations derived earlier via the BE identity. Let us 
comment on this difference. The lOj-symbol is not invariant under any Pachner move, thus it does not correspond to 
the quantization of a classical topological invariant theory. Furthermore, it is also not related to general relativity [35| . 
It does however capture some of its properties under the special restrictions of considering a single 4-simplex, and 
restricting attention to area observables [26l . [36l . [37j . Therefore the recurrence relation is implementing the symmetries 
of such toy model. 

However, this nice geometric interpretation of the recurrence relation comes from focusing on the Regge part of the 
asymptotics of the lOj-symbol and neglecting the other non-geometric terms [2^. The latter turn out to dominate the 
asymptotics: when the spins are scaled with a large parameter A, the Regge part goes as A~2, while the dominating 
constribution goes as and cor resp onds to degenerate configurations [2J, l2a]. Further intermediate non-geometric 
behaviours are also observed [13, [3J|. The situation is thus much richer than for the 6j: there, the geometric and 
non-geometric saddle points scale in the same way, and the second order recurrence relation Q includes them both. 
Here on the other hand the different scalings mean that recurrence relations will in general be of higher order. To see 
this, let us consider a simple example. We take the following asymptotic behaviour, 

{10j}~ ^ + j^ibcosXS + b'smXS) + ..., (48) 

where S is the (scale- free) Regge action for the 4-simplex (which defines the oscillation frequency of the lOj-symbol) , 
and a, b, b' are fixed numerical parameters. This form is motivated by the real asymptotics of the lOj symbol with 
the intermediate saddle points (see [HiIU) neglected. Then one can show that is a solution (for all parameters 
a, b, b') of a third order linear differential equation: 

A(/3(/3 - 1) - aS^\^){dl{\''+'^,f) + S^dx{\"+^f)) - ^{{{3 - l)(/3 - 2) - aS''X^){dl[\"+^ J) + S^{\"+^f)) = 0. (49) 

Extra terms in (|48p between the dominant term A^" and the oscillating regime will require a higher order equation. 
This example partially justifies the high order (fourth) of the recurrence relation for the lOj found above. However, 
we have not been able to study the asymptotics of our recurrence relation and make such a link with the saddle points 
of the lOj symbol more explicit. 

On the other hand, one might wonder whether restricting one's attention to geometric saddle points, second order 
recurrence relations with simpler geometric interpretations can be selected. There is indeed a strate gy t hat can be used 
to achieve this, and which has already proved fruitful in the computations of spin foam correlations [26l. [27l. [28l. [29l. [sst . 
It consists in adding a boundary state to the amplitude which naturally enhances the geometric saddle points, for 
instance suppressing the degenerate configurations in the bulk as in (26j . We develop this idea in the last section. 
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C. Geometric moves for the BC vertex 

Here we propose an interpretation of the different terms entering the relation above. The Barrett-Crane model 
is not topological, so that a geometric interpretation using Pachner moves, like for the 6js and 15js, is not possible. 
But, since the spins labelling the lOj-symbol are quantum areas, we would like to find some elementary moves which 
generate the area shifts appearing in the above recurrence relation. To this aim, the picture we obtained from studying 
the 2-4 move for 15js seems appropriate, because we interpreted there the move through the gluing of some flattened 

4- simplex, or the small displacement of a point of the initial 4-simplex, which geometrically results in some area shifts. 
We will show here that if we discard what happens at intertwiners in the recurrence formulas for 15j-symbols and 

focus on the behaviour of spins, then the different terms of the recurrence on the lOj look very much like terms which 
would arise if we were working on SU(2) 15j-symbols. It provides with a qualitative picture of the recurrence as the 
statement of invariance under a combination of the elementary move described in Section IIVBI in terms of gluing 
some flattened simplices to the initial one. However, in contrast with the 15j case, we will not give a precise algebraic 
translation of these geometric considerations. In particular, wc do not know the amplitude which should be assigned 
to the flattened simplices. 

We described two elementary moves in Section llVBI First, when slightly moving a point of a 4-simplex, the spins 
of three triangles sharing an edge are shifted by i, and we call it move 1 (it comes from taking A = ^ in the formulas 
of llVBp . Second, if the point is moved a little further (that is by taking A = 1), it is possible to shift the spin of only 
one triangle, by an amount of 1. This is move 2. Also remember that the links of the lOj graph stands for triangles 
of the 4-simplex, while vertices of the graph are tetrahedra. 

Let us begin with the simplest kind of shifts, when only one spin is changed by 0, ±1. This happens as soon as 
there is a permutation of two elements in (|l(7)l . Notice that no other triangle sharing an edge with the shifted triangle 
has its area changed. We can thus consider that such terms are obtained from the move 2. To make it precise, we 
use the situation described in section Hill with the same notations, and in particular the figures O and [SI We consider 
that the initial 4-simplex is (abcde), and wc proceed to a small displacement of the point a to a' - for 15j-symbols, 
the spins labelling (aa'c) and (aa'e) are taken to be zero, and that of (aa'd) is 1. Then, we have shown that it induces 
some shifts on the area of the triangle (abd) (again we do not look at the algebraic level and forget the shifts on 
intertwiners). Clearly, this elementary move can be independently performed on triangles which do not meet at an 
edge. 

A similar interpretation can be given to the terms consisting in a cyclic permutation of three elements. This induces 
shifts of ±i on three links forming a triangle in the spin network graph, or equivalently on the areas of three triangles 
sharing an edge in the original 4-simplex. This is precisely the move 1. If in the above-mentioned flattened 4-simplex 
(aa'cde) the spin labelling (aa'd) is ^ ~- or say, if it is smaller than for the move 2, since we do not have an algebraic 
characterization of the corresponding amplitude - then the three triangles sharing the edge (ad) are modified as 
expected. 

Consider now the terms due to a cyclic permutation of four elements in (|40p . thus acting on four triangles, whose 
corresponding links form a cycle in the spin network graph. If we think of a diagonal joining two opposite vertices of 
the cycle and dividing it into two triangles, we are in position to apply twice the move 1 for each cycle of three links. 
As far as 15j-symbols are concerned, some shifts have constant signs which can be arbitrarily chosen. This suggests 
that it is possible to leave invariant the spin of the chosen diagonal link when performing the move twice, while shifting 
the four links of the cycle. Let us precisely see how the successive two moves act on the original 4-simplex. Consider 
the 4-simplex of figure [5] and perform the move as described there. One obtains the 4-simplex (a'bcde), where the 
spins of (a'bd), (a'cd) and (a'de) differ from those of (abd), (acd) and (ade) by i. Add now a point a" close to a' and 
consider the 4-simplex {a'a"cde). Putting the spins of the triangles (a'a"c) and {a'a"d) to zero, and that of {a'a"e) to 
^, the spins of {a"be), {a"ce) and {a"de) in the new 4-simplex {a"bcde) differ from those of {a'be), [a'ce) (which are 
the same as those of {abe) and (ace)) and that of (a'de) by i. Between the triangles {ade) and {a"de), the area has 
been affected twice, in opposite ways, so that it is unchanged. Finally it is easy to check that the affected triangles 
{abd), {acd), {ace) and {abe) share the point a, so that they form a cycle of four links in the spin network graph. 

As for a cyclic permutation of the five elements in ((40|) . the strategy is to apply thrice move 1. In the spin network 
graph, it affects five links forming a closed loop (a 5-cycle). It is always possible to also consider two additional links 
so as to form three 3-cycles, such that two of them meet at a point and both are made with two of the links of interest 
and one of the additional links, while the third 3-cycle consists in the two additional links and the last link closing the 

5- cycle. We can then independently perform the move 1 twice to shift the six links forming the two triangles meeting 
at a point. The third move acts on the third 3-cycle, in order to eliminate the shifts on the additional links while 
changing the spin of the fifth link of interest. 
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D. "Ward" identities for BC correlations 

We consider the framework of the "spin foam graviton propagator" calculations on a single 4-simplex [2^, [l^, [H, [sgI . 
|37| . The boundary state is defined by the boundary probability distribution of the ten representation labels ji, .., jio- 
Here, for the sake of simplicity, we make the two usual assumptions: 

• that the boundary state is factorisable, i.e. of the type ni=iV'i(ji)5 where the ten functions -0^ are a priori 
arbitrary; 

• and that the boundary state is actually fully symmetric and "equilateral" , i.e. all the ten functions are the 
same, ipi{j) = '>Pij), Vi, where ip is a priori still arbitrary. 

For our present purpose, we do not need to require that the boundary state defined by be physical. We are 
interested in the correlation functions between the area of two triangles in the Barrett-Crane 4-simplex. To this 
purpose, choosing two triangles Ai = (ajS) and A2 = (7^), we define the following expectation value: 

(Oi(jaJO20aJ)v. ^ ^ E^i(-?aJ02(.7aJ n ^(j-f) ilOj}' (50) 

jab a<b 

where Oi and O2 are the two studied observables while the normalisation M is defined through the exact same 
sum without the observable insertions. To keep things simpler, we assume that the two observables arc the same, 
Oi - O2 = O. 

Let us now insert the recurrence relation on the lOj-symbol into such correlations. Shifts in the lOj-symbol can be 
re-absorbed through the summation over representation labels by a shift of the boundary states tjj and observables O. 
We will not write the full equation induced by the recurrence relation on the correlation functions, but we will focus 
on the structure of each term of that equation. 

Terms with shifts on the representations will involve shifted observable insertion such as 0(j ± ^) and ratios between 
the shifted boundary state and the actual boundary state ± ■^)/'4'ij)- The standard ansatz is to look at quadratic 
observables and to consider phased Gaussian boundary states: 

0{j) = f 0{j ± i) = 0{j) ± J + i (51) 
^(j) = e^d.^e-^^^-^"^"^^" ^^^^^ = e^^e-We*"('--), (52) 

V'(i) 

where i? and jo are the fixed parameters of the boundary states: i9 correspond to the classical (exterior) dihedral angle 
while Jo defines the scale factor of the 4-simplex (its size). And a defines the width of the Gaussian state. This way, 
we derive an equation between the 2-point function, the 3-point function, up to the 5-point function, with different 
observable insertions corresponding to shifts in the original observable or boundary state. It's a type of Ward identity 
for Barrett-Crane correlations on the 4-simplex. However, it is a priori not straightforward to interpret it as a function 
describing the scaling of the correlation function with changes of the scale factor jo. 



VI. CONCLUSION 



In this paper, we have derived recurrence relations for 6j-symbols and 15j-symbols. For 6j-symbols, we used 
different methods: the standard Biedenharn-EUiott identity which correponds to the 2-3 Pachner move, and the 
action of holonomy operators on spin network functionals. This second method naturally generates a regularized 1-4 
Pachner move and gives a precise relation to the classical constraints of the theory. Then, in 4d BF theory, recurrence 
relations for 15j-symbols have been derived from a regularized 2-4 Pachner move. The invariance of the Ponzano- 
Regge and Ooguri models under these Pachner moves is due to their underlying topological invariance, and provides 
the recurrence relations with a nice geometric picture: the shifts on the spins, which are eigenvalues of lengths and 
fluxes operators in the canonical theories, are generated by adding a flattened simplex to an initial one, resulting in 
a small displacement of a point of this simplex. We hope that such a picture will survive in non-topological theories 
and help to describe the action of the Hamiltonian constraints and to analyze the behavior of spin foam amplitudes 
under coarse-graining. 

Interestingly, the recurrences we found can be applied to any SU(2) {Snj {-symbol, depending on the form of their 
cycles. Indeed, we found that the relations naturally shift the spins along one cycle of the graph. Moreover, the 
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form and the coefficients of the relations only depend on the number of links of the cycle. We found relations for 
cycles made of up to six links. Besides its theoritical interest, we think that this feature could be useful to speed up 
numerical computations. 

For the special isosceles 6j-symbols, a new recurrence relation has been obtained using an integral representation. 
We extracted its geometric content by also deriving it from the BE identity. Then, we applied the same method to the 
lOj-symbol which is the 4-simplex amplitude of the Barrett-Crane model. While the recurrence formulas for the 6j and 
15j-symbols were difference equations of order 2, we found there a difference equation of order 4 for the lOj-symbol. 
This aspect is certainly related to the following fact. The evaluation of the lOj-symbol in the asymptotics, like that 
of the 6j-symbol, using the saddle point approximation j23 . |25| exhibits different kinds of saddle points: a geometric 
one together with non-geometric (degenerate) contributions. Although these different terms have the same scaling 
for the 6j, this is not the case anymore for the lOj, resulting in an increasing of the order of the difference equation. 
To get rid of the non-geometric saddle points, we then looked for recurrence equations for the vertex amplitude 
alone put together with a geometric boundary state and derived some Ward-like identities for insertions of geometric 
observables. 

We have also shown that the cosine of the Regge action solves the recurrence equation in the asymptotics. But we 
could not yet show that the other terms corresponding to degenerate saddle points also solve our recurrence formula, 
since we do not know an explicit formula for them. Looking at the Regge part of the asymptotics, the recurrence 
relation takes the form of a constraint imposing the closure of the 4-simplex. While the BC model is not topological, 
we have also proposed a qualitative picture to describe the various terms entering the recurrence, with the help of the 
elementary geometric moves which have emerged from the analysis of the 15j-symbol. 

Although the BC model probably does not provide a correct quantization of gravity, it is certainly a useful toy 
model to try to extend the results about spin foam models beyond the topological case. To complete the analysis of 
the lOj, it would be of great interest to derive the recurrence we found from a coarse-graining process, similarly to 
what we have done with the 15j-symbol. We also hope to apply the same methods to the more complicated Engle- 
Pereira-Rovelli & Freidel-Krasnov models [ij], which have a better semi-classical behavior. However, the intertwiners 
are then not frozen anymore and have to be taken into account in the recurrence relation and in the coarse-graining 
process. A first step in that direction could be to rewrite and reinterpret the 15j relations in term of the coherent 
state representation [l^ [s^ . 

APPENDIX A: THE USUAL RECURRENCE FOR THE {6j}-SYMBOL 

6j-symbols admit a recurrence relation of the following simple form: 

JlEUl + 1) L\ 1 I H + {2.11 + 1){2[J2(J2 + 1)^2(^2 + 1)+MJ3 + iMh + 1) - JlUl + l)h{h + 1)] 



[j2(j2 + l) + ?3(?3 + l)-Jl(jl + l)][;2(/2 + l)+j3(.?3 + l)-Jl(jl + l)]} [j.; '-J + + E{j,) |^.^'^^ ^. 



h h h 
h .7 3 



= 0, 
(Al) 



where: 



m 

It directly comes from taking one boundary spin to 1 in the Bicdenharn-EUiott identity, and evaluating explicitly 
the 6j-symbols which involve the spin 1. In the large spin regime, it can be turned into a second order differential 
equation, whose analysis leads to the known asymptotics of the 6j-symbol 

APPENDIX B: SOME REMARKS ON THE CHEBYSHEV POLYNOMIALS 



U2,{x) =2^'\{[x-cosJ^^. (B2) 



23 



APPENDIX C: GENERATING FUNCTIONAL FOR THE BC VERTEX 



We can also study the generating functional for the Barrett-Crane vertex. Indeed, using its expression in term of 
the Chebyshev polynomials C/„, we introduce the following function of 10 arguments tab for 1 < a < 6 < 5: 



[ 

ji6N/2a<6 a<b 



This is the integral of a rational fraction with a constraint. It looks like the evaluation of a Fcynman diagram and 
should be computable with similar methods (contour integrals). Or else finding a differential equation satisfied by 
-BC[t] should lead to recurrence relations on the lOj-symbols. Moreover, this generating functional could be directly 
relevant to the computation of correlation functions following the framework on the spinfoam graviton [26l . [27| . 
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